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Abstract 

Quantum computation has revolutionary potential for speeding algorithms and for simulating 
quantum systems such as molecules. We report here a quantum computer design that performs 
universal quantum computation within a single non-degenerate ground state protected from de- 
cohering noise by an energy gap that we argue is system-size-independent. Closely analogous 
to a traditional electric circuit, it substantially changes the requirements for quantum computer 
construction, easing measurement, timing, and heating problems. Using the standard adiabatic 
condition, we present evidence that this design permits "quantum concurrent processing" distribut- 
ing a quantum computation among extra qubits to perform a quantum algorithm of N gates in 
an amount of time that scales with ^/N. One consequence of our work is a fixed gap version of 
adiabatic quantum computation, which several arguments hinted could be impossible. 
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The standard gate model [ij for realizing a quantum computer involves fashioning a 
collection of two state quantum systems - qubits - that can be individually manipulated 
and measured. By pulsing the Hamiltonian of each qubit in a precisely timed fashion, one 
subjects the collection of qubits to a desired succession of unitary gates. If the experimental 
arrangement is capable of applying a sufficiently rich set of gates, universal computation is 
possible in which the initial state of the qubits |V'(^o)) can be taken along a time- dependent 
trajectory \ipiti)) , \4'it2)) , . . . to an arbitrary final state \ip{tN)). 

An alternate model, called ground-state quantum computation (GSQC) [2|-l5|, replaces 
the pulsed, time-dependent state with a "history" state |\E') in a larger Hilbert space 31 |. 
This history state possesses a time-independent record of the entire evolution of the time- 
dependent state, so that |\1') contains all of the amplitudes in \ip{tQ)), \ipiti)) , . . . , \ip{tj^)). 

GSQC works by establishing a Hamiltonian whose ground state is the history state. Let 
us start with the case of a single qubit possessing two quantum states |0) and |1) (see 
Fig. [D^A)). Suppose that a quantum algorithm specifies that our qubit begin at stage i = 
of the computation in the initial state |0) and undergo single qubit gates, given by the 
2x2 unitary operators Ui, i = 1, N (see Fig. lU^B) for A^ = 3). The state at stage i is thus 
supposed to be \ipiti)) = UiUi^i . . .Ui |0). We wish to capture in a ground state the history 
of this single qubit through all A^+1 computational stages. So, we define a Hilbert space with 
2(A^ + 1) states labelled by \bs), where 6 = 0, 1 is the bit value and s = 0, A^ is the compu- 



tational stage (see Fig. [HD)). Our history state is |\1') = J2i=o ivf^i • • • f^i |0i) /viV+T, a 
superposition of the system state at all A^ + 1 stages of the computation. It is useful to define 
in the 2(A^ + 1) dimensional Hilbert space a unitary operator Ui by {b'g>\Ui \bs) = {Ui)b',bSs,s'- 
Note that Ui acts on the bit value b of \bs) but does not change the stage variable s. To 
advance the stage variable, we define operators Agi^s = J2b=oi 1^*') that do not change 
the bit value b. The identity operator is /. Our history state can then be written |\E') = 



{I + An,n^iUn){I + An-i,n-2Un-i)---{I + AifiUi) |Oo) /VN+l. The state |^') is the ground 
state ofthe Hamiltonian H/e = |lo) {lo\ + J2b si^ ~^s,s-iUs) \^s-i) {bs-i\ (-^~^l^I,s-i) where 
e is the fixed energy scale of the Hamiltonian. This Hamiltonian is depicted schematically 
in Fig. [11(D) where a solid black line represents the on-site perturbation e |lo) (lo| and the 
dotted lines represent tunneling energies. There is an energy barrier to the excited states 
that scales like e/(A^ + 1)^ ^. 

Compared to the standard approach, GSQC requires a larger Hilbert space but demands 
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FIG. 1: (A) Schematic of Hilbert space of a qubit. (B) Standard gate model depiction of a single 
qubit undergoing an = 3 step calculation. (C) GSQC Hilbert space of a single qubit undergoing 
an A = 3 step calculation. (D) Hamiltonian for calculation in (C). Solid black line represents one- 
body onsite perturbation and dashed lines represent one-body tunneling energies. (E) Standard 
gate model depiction of quantum teleportation with simplified representation. (F) GSQC version 
of simplified representation of quantum teleportation. Solid wavy line represents two-body onsite 
interaction and dashed wavy line represents two-body tunneling. 



much less ti,n.depende„t control of a quantum system and provdes an energy barrier to 

decoherence [2] that defends against noise without active intervention. This is a fundamen- 
tal improvement over active quantum error correction schemes that constantly manage the 
qubits, detecting and correcting errors as they arise or echoing the qubits to average errors 
away. To see the value of this approach, note that with the exception of actual spin-1/2 
particles, most proposed realizations of quantum computers could not function without a 
protective energy gap separating two low energy qubit states from higher energy leakage 
states. GSQC takes this well-tested idea further, performing the entire multiple qubit com- 
putation within a single, non-degenerate ground state that is protected from noise by an 
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energy gap. In addition, the time-dependent pulsing Hamiltonian is discarded, allowing for 
less coupling to the outside world and less difficulty with heating and timing. Although de- 
parting from the standard gate model may seem unfamiliar initially, the GSQC history state 
is simply analogous to an ordinary electric circuit [4], with the input and output connected 
by a time-independent spatial network of wires. The standard gate model quantum com- 
puter is actually more exotic, behaving like an abacus, with qubits shifting back-and-forth 
in time like abacus beads. In other words, the standard gate model is like a mathematical 
calculation in which each the input to each step is overwritten by the output. GSQC is like 
a mathematical calculation in which all steps are shown on the page, so that each step can 
be checked against the next. We show here that not only can incorrect steps be penalized 
with an energy barrier but also that the history state makes it possible to perform some 
steps concurrently. 

We need to change the GSQC formalism as presented so far, since it suffers from at least 
two deficiencies if we wish to perform calculations with a large number of stages N. First, 
if is measured, the probability that it will collapse to stage s = N, where the answer 
UnUn-i---Ui \0n) is stored, decreases with like 1/{N + 1). Second, the energy gap from 
the computationally meaningful ground state to the ffist excited state decreases with A^, so 
that the energy barrier to decoherence becomes weaker and weaker. 

To solve both problems simultaneously, we observe that, from the standpoint of the final 
stage subspace {\0n) , IIa^)}, the probability N/{N + 1) that the history state is inhabiting 
another stage can be regarded as a sort of leakage. It is known that quantum teleportation [a, 
Q], shown in Fig. H^E), can be an effective means of combating leakage js]. One can perform 
teleportation in a GSQC context [^]. This does involve a complication since a quantum 
computer operator who performs a measurement in GSQC will collapse the history state 32| : 
the Bell measurement step of teleportation therefore requires non-trivial adaptation. Rather 
than simulating a Bell measurement with 4 possible outcomes, GSQC instead projects on 
to a singlet with only one possible outcome. This works as follows (see Fig. [H^F)). Imagine 
we give our initial qubit a 4 state Hilbert space \bs) where 6, s = 0, 1 and prepare it in the 
state (/ + Ai qUi) \0q) /v^, just as in the single qubit ground state quantum computation 
described above. Two qubits, each given a 2 state Hilbert space |6o), ^ = 0, 1, are in a Bell 
pair (|Oo) |lo) — |lo) |Oo))/v^- The three qubit state in the 2 (8> 2 (g) 4 Hilbert space is thus 
^(|Oo) |lo) — |lo) |Oo)):^(-^ + ^i.o^i) |Oo)- It is the ground state of the static Hamiltonian 
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'H + T-Lcreate pairs ) 

7^ = e/ ® 7 (g) |lo) (lo| + / / ® H{Ui), 

T^Create pairs p^ij-s ® 

H{U)/e = (1 - Ai,oW) \ho) (6o| (1 - U^A{,), (1) 

6 

Hcreate pairs 

A = 2 [(1°°^ 1°°^ + 1^°^ |lo))((0o| (Ool + (lol (lol) 
+(|lo) |Oo) + |Oo)|lo))((lo|(Oo| + (Oo|(lo|) 
+ (|0o) |Oo) - |lo) |lo))((0o| (Ool - (lol (lol)] ■ (2) 

We now incorporate the Bell projection which completes teleportation. We will need to 
store the post-projection states of the projected qubits and therefore supplement the Hilbert 
space, supplying a 5th state |IDLE) to the initial qubit and a 3rd state |IDLE) to the right 
member of the Bell pair. Thus, the net Hilbert space has dimension 2 (g) 3 (g) 5 instead of 
2 (g) 2 (g) 4. The (unnormahzed) ground state, including the Bell state projection, is 

I*) = i=(|0o) |lo) - |lo) |Oo))-^(/ + A,oWi) |0o) + |0o) IIDLE) |IDLE) . (3) 

Here, A is a real parameter that will be made large to increase the relative contribution of 
the post-teleportation state to |^). The Hamiltonian of the system is H + Hcreate Beii where 
T-Lcreate Bell IS Unchanged but now = e/ (g) / (g) |lo) (lo| + (g) H{Ui, A) where 

H(U, A)/e = 7 ® H{U)/e + |IDLE) (IDLE| ® l^s) {bs\ + Yl l^o) ® \IBLE) (IDLE| (4) 

6,5=0,1 6=0,1 

(^JO.) Hi) HI.) 10.) _ ijoLE) IIDLE)) (^A<~WA Jh^l _ (jo^EI {IDLE|) . 

Note that none of the matrix elements of Hamiltonian H{U,A)/e exceeds 1 irrespective of 
the value of A: there are no unphysical demands being made of our GSQC Hamiltonian. 

Assessing our GSQC teleportation procedure, we note that originally a GSQC of a single 
gate Ui acting on a single qubit produced a state (/ + Ai qUi) \0q) /a/2. Measurement 
extracted the answer Ui |0i) at the final stage s = 1 with probability 1/2 and the result 
|0o) at the earlier stage s = with probability 1/2. Measurement of the leftmost qubit in 
the new 2 (g) 3 (g) 5 dimensional Hilbert space now yields an answer at the final stage with 
certainty since this qubit has only 1 stage, but this answer is incorrect with probability 
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p = 6/(8 + A^). We have mapped the GSQC leakage problem into the usual standard 
model problem of faulty gates. Fortunately, this problem is well studied 9l-|l3| and one can 
eli,„,.ate the errors by applying a GSQC vers.cr of quantum error detecfng codesH- For 

instance, the GSQC version of the [[4, 1, 2]] Bacon-Shor code [8[ is described in the Appendix. 
Using the Hellman-Feynman theorem, one can argue that the energy gap to the first excited 
state satisfies the lower bound 

E{A) e-2arcsmh(A) ^ 2 



E{A = 0)- 1 + A2 4(1 + A2)2 Vl2. 

where E{A = 0) = (3 — \/5)e/2 (see Appendix). This lower bound should apply to an 
arbitrary number of GSQC qubits undergoing an arbitrary number of GSQC gates. In 
fact, this lower bound may underestimate the real gap substantially; small simulations even 
suggest E{A)/E{A = 0) ^ p rather than ^ [S^l- Note that the lower bound increases 
with error probability - a small error probability per gate requires a small gap. For a GSQC 
in the laboratory, the minimum achievable distortion p of the ground state is constrained 
by E{p) ^ ksT where fcsT is the temperature in units of energy. Thus, the usual quantum 
computing problem of reducing the gate error probability p below a threshold ^threshold has 
been exchanged for the potentially easier problem of having a sufficiently low temperature 
and/or high Hamiltonian energy scale e to allow -E(pthreshoid) ^ ksT. Note that, in addition 
to the inherent error probability 6/(8 + A^), imperfections in fabrication must be added to 
p when determining if p < ^threshold- Such time- independent errors can be minimized by 
testing and refining the apparatus and by compensating for known distortions using GSQC 
versions of decoherence free subspaces 14, 15|, dynamical decoupling |l6|, and the like, while 
both time-dependent and time-independent errors will be suppressed by the energy gap (see 
Appendix) . 

It is straightforward to generalize to a single qubit calculation with operators 
Ui, . . . ,Un, as shown in Fig. [2] for the case A^ = 3. The Hilbert space is a tensor prod- 
uct 2(8)3(S)5(8)---(8)3(S)5(8)3---(8>5. The generalization to multiple qubit calculations is 
similarly straightforward (Appendix). This GSQC design requires an amount of hardware 
that scales linearly in the number of gates, like a traditional electric circuit with no clock 
cycle. 

All of the steps of the calculation are simultaneously present in the GSQC history state, 
and GSQC appears to be capable of a new type of quantum speedup we term "quantum 
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FIG. 2: Iterated GSQC quantum teleportation that performs = 3 step calculation of Fig. [HD). 



concurrent processing." Our argument employs adiabatic quantum computation (AQC). 
AQC and GSQC both involve quantum computation using ground states, but the theories 



have different origins. In 2000, Farhi et a. 



suggested using adiabatic evolution to solve 



"instances of the satisfiability problem" 17|, |l8|. They continued to develop AQC as a cre- 
ative methodology for attacking certain types of computational problems but not a form 
of universal computation that could execute existing algorithms. In contrast, when it was 
proposed in 1999, GSQC had no new algorithmic content. Instead, GSQC provided a uni- 
versal history state and Hamiltonian suitable for computing, for instance, Shor's algorithm 
[lol . while reducing the need for time-dependent control and providing an energy gap for 
combating decoherence. AQC gave an adiabatic recipe for reaching the ground state which 
they carefully distinguished from cooling 17|]; the adiabatic recipe can be justified using the 
adiabatic theorem. GSQC, on the other hand, did imagine cooling to the ground state but 
never analyzed whether this would work for the GSQC Hamiltonian 35|]. After AQC was 
proposed, Mitchell realized it was probably straightforward to reach a GSQC ground state 



by adiabatically tuning the parameters of the GSQC Hamiltonian 



36|, albeit not necessar- 
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ily in the canonical Hq{1 — s) + His fashion of AQC, and an exphcit formulation of this 
"squeegee" method has appeared [sl. And since GSQC was proposed, AQC researchers have 
investigated the efficacy of an energy barrier to decoherence j20l] and have come to embrace 
the idea of a universal history state [2l| . 

A number of recent developments have given the impression that the AQC energy gap 



shrinks quickly with system size, leading to pessimism |22|-|24| about AQC since a small gap 
translates into unacceptably long AQC calculation time and vulnerability to environmental 
excitations. The version of GSQC formulated in this paper, however, yields a fixed-gap 
version of AQC as a corollary 37|. To demonstrate fixed-gap AQC, one need only show that 
it is possible to tune the parameters of the fixed-gap GSQC Hamiltonian adiabatically to 
bring the system into the GSQC ground state. This can be achieved by starting with A = 
in the GSQC Hamiltonian. The GSQC Hamiltonian is then extremely simple, essentially 
organizing the qubits into Bell pairs. The gap is (3 — -\/5)e/2. Next, we adiabatically take 
A from to A^az at a rate that depends on the system size. The traditional adiabatic 
condition states that the system will remain in the ground state as long as the increase of 
A occurs on a time 

T » 1(^,1 dU/dK m /El (6) 

where k is any excited eigenstate of the net Hamiltonian and we have simplified the denom- 
inator using the fact that the energy of vanishes. As is shown in the Appendix, we are 
guaranteed to satisfy this inequality if T ^ hKmax \/ Ne / E'^ ( A^ax ) , where E^Arnax) is given 



38 



by ([5]) and is the number of computational steps 

Refiecting on this finding, we see that it implies that GSQC allows the computation of 
an N step quantum algorithm in 0{\/N) steps (not including a 0(poly(log(A^))) factor for 
error correction). For instance, consider a computation consisting of a series of black-box 
single qubit gates. Since the gates are black boxes, we can devise no clever shortcuts to 
speed the algorithm; in the standard gate model, we would need 0{N) time to apply the 
gates in series. Using GSQC, the result of this series can be obtained in a time 0{\/N) 



using 0{N) resources 391]. Since the calculation is distributed among 0{N) times as much 
hardware, we term this result "quantum concurrent processing." It is quite distinct from 
the familiar "quantum parallelism" that requires clever algorithmic design; here there is no 
structure at all assumed about the algorithm. It is also distinct from Grover's search of 
an unsorted database of N items with 0{\fN) oracle calls - there is no oracle involved in 



quantum concurrent processing. Evidently, a degree of coordination of resources is possible 
in GSQC that is not possible classically; using entanglement, there is no need to wait for 
the first part of a task to be completed before the next part of the task gets started. 

We gratefully acknowledge helpful comments from Marvin Kruger, Daniel Lidar, Keith 
Miller, Kevin Osborn, and Mark Wilde. 
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Appendix 

Universal quantum computation 

We have used first quantized notation in the main text to render our findings as accessible 
as possible, but this notation becomes unwieldy as the number of gates and qubits increases. 
Therefore, we now return to the second quantized notation of Refs. j2"5|. In this notation, 
the state \bs) of qubit Q is expressed using a creation operator Cg^^lvac) where Q labels 
the qubit, s labels the stage, and b labels the bit value. Grouping the creation operators 



into row vectors, we define Cq 



. For simplicity, we assume that the qubit 



particles have bosonic commutation relations, although this does not affect any results. In 
this notation, the stage operator Ag/^g = J2b=oi 1^*') for qubit Q just becomes Aqs'^s = 

Substituting the second quantized notation into the history state = (/ + 



^iV,7V-lW^)(/+^iV-l,^-2W,V-l)...(/+^l,oWi) |Oo) /VNTT = E^=0,...,7V ■ ■ ■ Ul \0i) /VnTT 

for qubit Q yields = (/ + Cj^^t/Q,7vCQ,7v-i)(/ + Cj,;v-i^Q,JV-iC^Q,iv-2)...(/ + 

1 







|vac) / \/N + 1. The Hamiltonian of this qubit T-L/e = \1q) (1o| + 



Zbsi^-^s,s-iUs) \bs-i) {bs-i\ {I-UlAl,_,) can be written as H'^/e = C^^ 




1 



We now describe the GSQC teleportation process that applies Uq^i to qubit Q. First, 
consider the case of a single gate f/g,!, the case which was described in the main text and 
depicted in Fig. lU^F). Label the three particles involved in the teleportation Qi,Q2, and Q3. 
Set the |IDLE) states of Qi and Q2 to rfg^ |vac) and rfg^ |vac) respectively. The Hamiltonian 
is then "H + V-create Bell whcrc 



n = eC, 



Qi,0 




1 



T~i-Create pairs -^Create pairs 



(7) 
(8) 
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i^(Qi,Q2,f/Q,i, A)/e = - CI^,,oUI,){Cq,,i - f/Q,iC^Q„o) 

and 



Q2"Qi 



(Cq3,0,oCq2,0,0 + Cq3,0,i'^Q2,0,i)('^Q2,0,0CQ3,0,0 + CQ2,o,iCQ3,o,i) 



(cq3,o,iCq2,o,o + '^Q3,o,oCq2,o,i)(cQ2,o,oCq3,o,i + CQ2,o,iCQ3,0,o) 
+ (Cq3,0,0^Q2,0,0 ~ c^3,o,iCq2,o,i)('^Q2'0'O^Q3,o,o - CQ2,0,iCQ3,0,i) 

The ground state of this Hamiltonian is 

'^Q3.o.o'^Q2,o,i ~ '^Q3.o.i'^Q2,o,o ^t 



(10) 



I vac) 



/^(-^ + ^Qi,l^QA'^Qifi)^Qi,0 



+ 



A 



2V2 



Cq.^,oUq,i 



I vac) . 



With the second quantized notation, it is straightforward to generahze to treat the case 
of a single qubit Q undergoing a succession of gates Uq^i, i = I, ...,N (shown in Fig. |2]for 
the case = 3). It has Hamihonian T-L + Ticreate Beii where 




1 



N-l 



CQi,0 + ^ H{Q2q+l, Q2q+2, Uq^g+l, A) 
q=0 



and 



N-l 



'Hcreate pairs 

q=0 



(11) 



(12) 
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in terms of (|9]) and (fTOj) . Its ground state has the form 

tN- 



1^) = nf-^ 



(I + A(i"^ '^Q2q+2,0fiCQ2q + lA,l - CQ2q+2,0,lCQ2q+l,l,0 -. 

^ Q2q + 2 Q2q + 1 

+ '^Q2,+i,i^Q,g+iCQ2,+i,o) 



TiN-l ^Q2q+3,0,0^Q2q+2,0,l '^Q2g+3,0,l'^Q29+2,0,0 



I vac) 



(13) 



So far, we have addressed a single qubit Q undergoing a series of single qubit gates 
Ug^i with teleportation steps in between the gates. Of course, for two qubits Q and R each 
independently undergoing a series of single qubit gates f/g^j and Ujij, the ground state would 
take the product form 



l*> = n,".-,' 



n 



7V-1 
r=0 



Q2(j+2 <92g+l 

,t C/J2^^2_0,0C_R2r+l,l,l ~ <^^?2r+2,0,lC_R2,.+ i,l,0N 1 



(/ + Adi di 



72 



72 



^Ar_l '^Q2q+3,0,0'^Q2g+2,0,l ^Q2<;+3,0,l'^Q29+2,0,0 ^t 



^Af-l '^kr + 3.0.0'^kr+2,0,l + 3 .0. 1 '^Q2r + 2 .0,0 

r=0 /TT ^-Ri.O 



72 



1 



1 





I vac) 



(14) 



and the Hamiltonian would simply be the sum of 

N-l 




1 



+^^k,o 




1 



CQi,0 + ^ H{Q2q+l, Q2q+2, f^Q,g+l> ^) 
9=0 

N-1 

Cr^A + H{R2r+l, R2r+2, f^i?,r+l5 A) 



r=0 



and 



TV-l AT-l 
'Hcreate pairs = Hcreate Beu{Q2g+2, Q2q+z) + Hcreate Bell{R2r+2-, R2t+^)- 



q=0 r=0 

A non-trivial quantum computation will require two qubit gates. To include a controUed- 
PHASE gate between qubits Q and R at stage s, we remove the g = s — 1, r = s — 1 factor 
in the product state f|T4l) . which is 
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We replace this factor with 



(I + Ad) '^Q2s,o,o'^Q2.-i,i,i ~ Cq3^,o,i'^Q2.-i,i,o ^ 

^ Q2s Q2S-I 

(I + Ad^ £R2^£^oCij2^_i,i,i - c/j2^_o,iCfi2,_i,i,o^ 



1 

V2 



~7k(^ + cks-l,l,0Ci?2.-l,0,0CQ2,_i,lCQ2.-l,0 + ck,_i,l,lC«2s-l,0,lCQ2,_i,l ^C'q2s-1,0) 



where Z is the 2x2 Pauh matrix o"^. The quantum state is now no longer a product of the 
state of qubit Q and the state of qubit R since the controUed-PHASE gate entangles them. 

To make this the ground state, we must alter the q = s — 1 and r = s — 1 terms 
of the Hamiltonian, H{Q2s-i,Q2s,Uq^s, A) + H{R2s-i, R2s,Ufi^s, A). We replace the single 
qubit gate Hamiltonians (C^^.-i,! " C'q2s-i,o^q,J(<^Q2s-i,i - Uq^sCq^^.^o) and (C^^.,,! - 
Ck.-i,of^L)(*^«2.-i,i - f^i?,sCij2,_i,o) with a two qubit gate 

(cks-l,l,l*^Q2.-l,l ~ ck^-l,0,l*^Q2s-l,0^)('^Q2s-l,lCi?2s-l,l,l - ^C'Q2s-l,0Ci?2s-l,0,l) 
+ (•^1223-1,1,0^023-1,1 ~ '^k3-l,0,O^Q23-l,o)(^Q23-l,lCiJ2s-l,l,0 - C'Q2s-l,0CiJ2.-l,0,0) 
+ (^ks-l^«2.-l + C'ks-l,l*^«2s-l,l)CQ2s-l,0'^Q2s-l,0 

+(^Q2s-i^Q2s-i + C'q2s-i, 1*^02^-1, i)Cks-i,oC'i?2.-i,o- (15) 

The last line two lines impose an energy penalty if one of the qubits crosses through the 
gate without the other. The CPHASE gate is depicted in Fig. [31 

Since we have considered both single qubit and CPHASE gates, we have a universal 



construction [25| for universal GSQC with teleportation. Naturally, we could use a CNOT 
instead of a CPHASE by replacing the 2 x 2 Z matrices in ( ITSll with 2 x 2 X matrices. More 
generally, we can write 

(CL-1,1 ® C'^2.-l,l - 42.-1,0 ® C^L-l,0f^^)(C/?2.-l,l ® C'o.^-l,! - ^C'fi2.-1,0 ® C'Q2.-1,o) 
+ ('^ks-l^'R2s-l + C'ks-l, 1*^^2.-1, i)C'q2.-1,0'^Q2s-1,0 
+ ('^Q2s-l^Q2s-l + C'Q2s-l,l'^Q2s-l,l)Cks-l,0C'H2s-l,0 

where is a 4 x 4 two qubit unitary matrix. (In the CPHASE case, for instance, U = 
f^cPHASE = diag(l, 1, 1, -1).) 

A two qubit Grover [2(y\ algorithm provides an illustrative example of a GSQC compu- 
tation. This algorithm has circuit depth = 5: the two qubits start in the product state 

13 



^25+1 



o 

o> 

OStJ 



^2s -..( )IIDLE> 



5?- 



O^G .niiDLE> 



O ^ 

O 
V) 



.2s - - - ''IDLE> 



>;2s+i 



FIG. 3: Schematic depiction of CPHASE gate. 



|0) |0), undergo Hadamard gates H, interact via an oracle (one of the 4 possible controlled 
phase gates), undergo Hadamard gates, interact via the phase gate I ® I — 2 |00) (00 1, and 
undergo Hadamard gates. The result is the computational basis state whose phase was 
flipped by the oracle. The algorithm is shown in Fig. IHin standard gate model, GSQC, and 
teleportation-based GSQC versions. 
Energy Gap 

To compute a lower bound for the energy gap, we employ the Hellman-Feynman theo- 
rem. This theorem states that for a parameter dependent Hamiltonian H{A) with a non- 
degenerate eigenstate |\E'fc(A)) whose energy is -E'fc(A), 

^ = (^^(A)I^I^.(A)). (16) 
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FIG. 4: (A) Two qubit Grover algorithm in the standard gate model. (B) Two qubit Grover 
algorithm in GSQC without teleportation. This version of GSQC is not scalable since the gap 
decreases with circuit depth N. (C) Two qubit Grover algorithm in fixed gap, teleportation version 
of GSQC. Each qubit has a snaking geometry, as depicted in the small flow diagram. Each gate 
requires a Hilbert space of size 2(8)3(X'5 (e.g. 3 electrons in 10 dots). The fabrication requirements 
are thus linear in the number of gates, just as in an ordinary electrical circuit. 



Now, differentiating ([9]) gives 
2 



^ CQ„o,oCQ„l,l - CQ,o,lCQ„1.0 _ ^^^^^^ , 1^^^ 



V2 



1 + A2 



•^Qi, 1,1*^(32,0,0 '^Qi,l,0'^Q2,0,l 



V2 



+ HA2 dQ2dQi 



1 



1 + A2 
2 

~1 + A2 



^A^ '^Qi.l. 1*^(32,0,0 '^Qi,l,o'^Q2,0,l _ ] ^^A^ '^(32.0,0(^Qi,l,l - CQ2,0,lCQi,ifl _ 



V2 



^Qi, 1,1*^02,0,0 '^Qi, 1,0*^02,0,1 



V2 



V2 



+ ^^01^2 dQ^dQ^{Ek-h{Q,,Q2))\'^k) 
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where h{Qi,Q2) is defined as the sum of terms in the Hamiltonian that do not involve Qi 
or Q2 and where we have used the Schrodinger equation for |\I/fc). Now, 



1 + A2 



t _ 



•^gi,!, 1*^(32,0,0 '-Qi,l,0'-^Q2,0,l 



+"Qi"Q2 ^ ^^"Q2"Qi 

since the smallest eigenvalue of the matrix — [ ? 2A ] is —A — y/l + Is?. Note that we have not 
proven this last inequality in all cases; this is why we use the term " argue" rather than the 
term "prove" in the main text and in the abstract. Our final result is obtained by summing 
over Qi and Qj 



9=0 

where the 2 in the final equation arises because interaction terms between Q2q+2 and Q2q+2, 
appear in neither h{Q2q+i,Q2q+2) nor /i(Q2<?+3, Q2g+4)- 

This differential equation for Ek{A) has the solution given in the main text, 

E{A) e-2arcsmh(A) 

E{A = 0) - 1 + A2 ■ ^^^^ 

When A = 0, it is straightforward to obtain the eigenvalues of the Hamiltonian exactly. 
Error correcting codes in GSQC 

We established in the main text that GSQC allows fixed gap computation provided one 
allows a non-zero error probability per gate. Fortunately, the problem of handling faulty 
gates is well studied in the context of the standard gate model; we need only adapt the con- 
clusions of fault tolerant quantum error correction to GSQC. Some complications do arise 
from the properties of measurement in GSQC. One can only perform a measurement at the 
end of a GSQC, not in the middle of a calculation. It is known that quantum error cor- 
rection can be performed coherently without measurement, and this could be implemented 
in GSQC. However, the noise threshold is expected to to decrease markedly if one refrains 
from measuring j27|. 
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Fortunately, as we showed in the main text in the case of Bell measurements, in GSQC 
one can modify the Hamiltonian to project j4| a single result in the wavefunction rather 
than performing a genuine measurement with different possible outcomes. This is partic- 
ularly powerful in the case of error correction, since one can project the error syndrome 
measurement to a "no error" result. (This in turn makes it possible to use error detecting 
codes rather than error correcting codes for GSQC if desired, since one always projects to a 
"no error" result and does not need to apply corrections.) 

To see how this projection works, suppose that our data qubits have GSQC ground state 
at some stage s of the calculation. We interact them with an ancilla to obtain the 
error syndrome, and, after the interaction, the zero-energy ground state of the Hamiltonian 
is |Oo) l^&o) + |lo) l^i)- Here, the ancilla state |Oo) indicates a "no error" syndrome, and 
|\l/o) is the error-free state of the data qubits, while |lo) indicates a non-trivial error in 
the state of the data qubits. In the standard gate model, we would measure the 
ancilla and apply a correction to the data if the result |lo) were obtained. In GSQC, we 
instead amplify the |Oo) part of the wavefunction and never apply a correction. To achieve 
this, we supplement the Hilbert space of the ancilla qubit with an |IDLE) state and add 
a term to the Hamiltonian of the form e(A |Oo) - |IDLE))(A (Oo| - (IDLE|)/(1 + A^) or 
e(AcQQQ — dQ){AcQfi^Q — dQ)/{l + A^) in second quantized notation. Then, the new zero- 
energy ground state is [(|Oo) + A |IDLE)) |^o) + |lo) l^i)]/(l + A^). The "no error" result 
|\l/o) has now been amplified; it can be made very large by increasing A. 

We apply this projection method to a GSQC version of the smallest Bacon-Shor error 
detecting code. This is a distance 2 subsystem code that encodes one logical qubit into four 
physical qubits (i.e. it is a [[4,1,2]] code). Usually, the four physical qubits are imagined 
to lie in a square with coordinates (1,1), (1,2), (2,1), and (2,2), but we will label them 
Q, R, S, and T respectively. In the standard gate model, the stabilizer generators can 



be measured 



28 



via the circuits shown in Fig. |5l It is straightforward to translate the 
CNOT gates appearing in these circuits into GSQC; this was described above. We translate 
the Z measurements into Hamiltonian terms e(A |Oo) - |IDLE))(A (Oo| - (IDLE|)/(1 + A^) 
as described above and the X measurements into Hamiltonian terms e(A(|Oo) + |lo))/v^ — 
|IDLE))(A((Oo| + (lo|)/y2-(IDLE|)/(l + A2). This projection amplifies the no-error portion 
of the wavefunction, obviating the need for a correction circuit. A schematic of the GSQC 
version of the lower left circuit in Fig. O is shown in Fig. El To make a rough underestimate 
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FIG. 5: Error detection circuits for Bacon-Shor [[4,1,2]] code from 28(. The four physical qubits 
are labehed Q, R, S, and T. 



of the pseudo-threshold, we notice that there are 12 gates in the circuits in Fig. [5l including 
measurements since they are translated into imperfect projections. A CNOT extended 
rectangle will therefore have 4 x 12 error detection gates, and 4 transverse CNOT gates. 
This is a total of 52 gates. If 2 gates fail, the error is undetectable and cannot be projected 
out, and there are about (^^^j, or about 1000 ways for this to occur. The resulting pseudo- 
threshold estimate is 1/1000 = .1%. Note that (i) this pseudo-threshold is quite high 
compared to most in the standard gate model literature and could be an underestimate by 
an order of magnitude given that we didn't distinguish between malignant and benign faults 
{40! , (ii) the overhead involved in our error correction is exceedingly low, requiring very few 
extra qubits and no classical computation abilities in stark contrast to most high-threshold 
standard gate model error correction schemes, and (iii) this subsystem code is well-suited 
to locality constraints. We considered the [[4,1,2]] code for simplicity; it is very likely that 
qualitative improvements will be possible by considering other existing codes and/or refining 
existing codes to take advantage of the properties of GSQC. 
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FIG. 6: GSQC version of lower left error detection circuit in Fig[5l 



We emphasize that this error correction method will handle fabrication errors as well 
as the errors introduced by teleportation as long as the net probability of error is below 
threshold. On the other hand, time- dependent perturbations of the system are opposed by 
the energy gap without any active error-correction. A time-dependent perturbation of the 
form Hpert = Vpert COS ut cau introduce an oscillating contribution Cfc(t) |\E'fc) from the excited 
state into the ground state wavefunction. Time-dependent perturbation theory gives 



^i{EklTi+u)t _ ^i{Ek/h-u))t _ ^ 

+ 



Ek/h + LO Ek/h — LO 

If the excitation frequency is far below the gap frequency, this coefficient will not become 
substantial, and the measured GSQC wavefunction will not suffer distortions. For excitations 
with frequencies approaching the gap frequency, transitions out of the ground state can occur 
at a rate approximated by Fermi's golden rule 

Adiabatic calculation 
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The adiabatic condition states that system in the ground state of the Hamiltonian will 
remain in the ground state if all of the parameters in the Hamiltonian change sufficiently 
slowly. There have been some recent attempts to make rigorous adiabatic conditions, es- 
pecially in the context of adiabatic quantum computing, but equation in the main text 
gives the traditional condition for our problem: 



T:^hA^a.\{'^k\dn/dA\^)\/El 



(19) 



where T is the time to take A/A^ax from to 1, |\&) = |\1'(A)) is the ground state of the 
Hamiltonian, k labels any excited eigenstate of the Hamiltonian, and we have simplified the 
denominator using the fact that the energy of |\E') vanishes. To evaluate the right hand side, 
we can use the identity 

U -C/ground 

A; 7^ ground 



72z? . ^ l/iTf. I ^17 /^A liTr\|2 ^2 



dA2 



Eh 



dK^ 



Since -Eground = for all A, this implies that for any term in the sum 



{^!u\d'H/dAm^ im'^m 



Noting that 



dA^ 



El 



2 El 



2\3 



'^Qs,l,l'^Qs+i,0,0 '^Qs,l,0'^Q«+i,0,l 

''5s Qs+i 

2-6A2 6A-2A3 
6A-2A3 2-6A2 

Qs+i 



/(1 + A^) 



we find that (^'| ^ |^) < 2A^e/(l + A^). The resuh is that it is sufficient to have 



A, 



Ne 



E{K 



where E{Amax) is the minimum energy gap, calculated above in the energy gap section. 
Note that here is the total number of gates in the circuit. If there are many qubits 
undergoing gates independently, it may be advantageous to adopt a more nuanced schedule 
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for increasing the A of the gates from to hmax- For instance, it may be advantageous to 
keep some A fixed at zero temporarily or to use the strategy that appears in js]. 

One intriguing theoretical problem is to determine whether an adiabatic approach to the 
ground state is necessary or whether a cooling process is likely to find the ground state 
quickly. 
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